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Since H 2 S was detected in Comets Austin (1989c 1) and Levy 
(1990c) in the microwave range, there has been increasing interest 
in searching for SH, which is the prime dissociative product of 
H 2 S. We present g - factors for the A-X (0-0) band of SH as a 
function of heliocentric velocity at r = L0 AU. We derive an 
upper limit production rate, g(SH)/g(H 2 0) < 0.017, for Comet 
Brorsen-Metcalf (1989o) and calculate a dissociative lifetime of 
105 sec at a heliocentric distance, r = 1.0 AU, and at a heliocentric 
velocity, v(r) = -28.5 km sec" 1 . *'■ 1992 Academic Press, Inc. 


Introduction . Colom et al. ( 1990) recently reported the first detection 
of H 2 S in the millimeter-wavelength range of Comet Austin ( 1989c 1). 
Subsequently Crovisier et al. (1991) presented production rates of H 2 S 
for Comets Austin and Levy (1990c) derived from their microwave data. 
Marconi et al. (1990) first identified H^S* in a mass spectrum of Comet 
Halley obtained with the heavy ion analyzer (PICCA) on board the 
Giotto spacecraft and attributed its existence to protonation of H 2 S. It 
is. therefore, of interest to search for SH , which is the prime dissociative 
product of H 2 S. Krishna Swamy and Wallis (1986) and Wallis and 
Krishna Swamy ( 1987) tentatively identified the B-X bands of SH in the 
1670-1920 A range in the IUE spectra of faint comets and of Comet 
Halley, respectively. The suspected features in those spectra are, how- 
ever, hampered by the typical IUE noise (e.g., Kim and A’Hearn 1991 ). 
The recent availability of good, high-resolution, groundbased spectra 
covering the spectral range including the SH A-X transitions (O'Dell et 
al. 1990) provides additional stimulus for estimating the expected 
strength of the emission by SH. 

g-factors of the SH A-X {0-0) hand. The individual lines of the SH 
A-X (0-0) band have been observed in absorption, and their line posi- 
tions were presented by Ramsay (1952). We used his line positions 
in our ^-factor calculations because the line positions are sufficiently 
accurate for calculation of the Swings effect, and they are the best 
available in the literature. 

The A-X system is the same type of electronic transition as is the OH 
A-X system. For the Honl-London (H-L) factors of this band, we 
adopted an updated list of the H-L factors for the - : II system of 
the OH A-X band presented in Appendix B of Schleicher and A’Hearn 
(1982). 

Friedl et al. (1983) found experimentally that predissociation of SH 
occurs throughout the rotational levels of the ()-() band of the A-X 
system, although previously only the v' s 1 levels of the A state had 


been thought to predissociate significantly. They found that the predis- 
sociative and radiative lifetimes are 3 and 820 nsec, respectively. Sene- 
kowitsch et al. (1985) subsequently revised the radiative lifetime to be 
704 nsec. Since most of the radicals predissociate subsequent to their 
first excitation to the A state, fluorescent equilibrium will not be estab- 
lished and the emission spectrum w ill be controlled by the Swings effect 
operating on a ground-state population. The general formula for the 
^-factor (,v ') for a transition undergoing predissociation can be ex- 
pressed by 

Hff = yjr X -v« r B rrPn • < 1 > 

r 

where the summation is restricted by the selection rule for the rotational 
transitions. ,v or is the fractional population of the/" rotational level in 
the ground state, y j r is the branching ratio for radiative decay from the 
/ rotational level of the excited state to the /' rotational level of the 
ground state, B ri is the Einstein B coefficient, and p r is the solar 
radiation density. There is major uncertainty as to what should be used 
for the populations of the ground-state rotational levels, .Vq,~. The follow r - 
ing is a discussion of the possible population of the ground state and 
corresponding temperature of SH. 

Hawkins and Houston (1980) conducted laboratory experiments on 
the photodissociation of H : S and measured the rotational temperatures 
of SH. They found that the SH radicals are formed a little warmer 
(about 40 K greater) than their parents when their parents are al room 
temperature. Hawkins and Houston (1982) performed the same experi- 
ment for H : S al 3 K and found that the rotational temperature of the 
nascent SH is about 150 K. This indicates that if the parents are very 
cold, the nascent SH becomes much warmer (150 K) than the parents. 
From the results of Hawkins and Houston, the rotational temperature 
of the nascent SH seems to be greater than 150 K regardless of the H : S 
temperature. 

Crovisier et al. (1991 ) measured a rotational temperature of 40 K for 
H 2 S in Comet Austin at about I AU heliocentric distance. Although the 
rotational temperature was not measured, we expect a low rotational 
temperature for this molecule in Comet B-M because rolationally ex- 
cited H 2 S relaxes rapidly (D. Bockelee-Morvan, private communication, 
1991). From these arguments, we can say that the temperature of the 
nascent SH is greater than 150 K and is likely to be less than 400 K 
because of the probably low excitation temperature of H 2 S. 

In order to estimate the changes in the rotational populations of SH 
expected after formation but prior to fluorescence, we roughly calculated 
rotational transition probabilities of SH as follows. The OH and SH 
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ground states have the same transition types, and therefore the only 
major differences are the frequency and dipole moment. The radiative 
transition probability is proportional to the square of the dipole moment 
and the cube of the frequency. The dipole moments of SH and OH are 
0.7580 and 1.6676 D, respectively, and the rotational B constants are 
9.46 and 18.5 cm" 1 , respectively. Burdyuzha and Varshalovich (1973) 
presented Einstein A coefficients for the hyperfine transitions of OH. 
We scale down these values using the ratios of the dipole moments and 
the frequencies of SH and OH. The resultant radiative lifetimes for the 
lowest lines in the H branch are about 200 sec, which are greater than 
the SH lifetime of ~ 100 sec as discussed in the next section. The radiative 
lifetime decreases as rotational quantum number J increases, and it 
approaches 1 sec at around J = 25. The radiative lifetimes for lines in 
other satellite branches are much greater than the SH lifetime. We 
conclude that the lower rotational levels of the SH radicals do not 
undergo significant rotational deexcitation (or cooling) before they are 
destroyed by photodissociation. 

Since the SH is formed where densities in the coma are moderately 
high, collisions can also affect the distribution of populations. The recent 
model of Crifo (1991) predicts kinetic temperatures near 200 K at a 
distance of 4000 km from the nucleus, the mean distance at which H,S 
photodissociates to produce the SH. Thus collisions, which would tend 
to produce rotational distributions in equilibrium with the kinetic temper- 
ature, will tend to produce a rotational excitation characterized by a 
temperature near 200 K. We therefore will assume that the rotational 
distribution in the ground state can be characterized by temperatures 
between 150 and 400 K. This range leads to about a 15% variation in the 
^-factor for the entire 0-0 band and about a 50% variation for the R ] + 
r Q 2] branch lines clustered near 3238 A. These variations are still smaller 
than the variation caused by the Swings effect. The variation for individ- 
ual lines within the band, however, is greater than that caused by the 
Swings effect. 

The other parameters needed to solve Eq. (1) are much belter deter- 
mined. The branching ratio, y, T , can be expressed in terms of the Ein- 
stein A coefficients for the downward transitions coupled with the disso- 
ciative lifetime of the upper state as follows: 


f ' n 'vT 

J > 1 l/'T'j) + Ajy + Ajj) i -I" Ayy _ | 


Here Ayy is the Einstein A coefficient and T d is the dissociative lifetime 
of the u' - 0 state in the A state. The equation is written for a {7-branch 
line, and the subscripts in the denominator differ for P and R branches. 
In all cases, the branching ratio can be approximated as Ayyr^ because 
1/T d dominates the term in the denominator. 

% is the Einstein coefficient for radiative absorption in a given 
transition and can be expressed in terms of A(0-0), the Einstein coeffi- 
cient for spontaneous emission in the entire band, as follows: 


B = A(0-0) H-L 

JJ 8t7/?wJ_„(2/" + 1)' 


(3) 


where h is Planck's constant. oj 0 _ 0 is the frequency of the band head in 
units of cm 1 , and H-L is the Honl-London factor for the transition 
taken, as noted above, from Schleicher and A'Hearn (1982). 

The radiation density at each transition, p rr , is taken directly from 
A’Hearn et al . (1983) with allowance for the Swings effect. 

Using the above equations, ^-factors for Comet B-M (1989o) have 
been calculated for a heliocentric distance, r - 1 .0 AU, and for a 
heliocentric velocity, v(r) = -28.5 km sec -1 . A model spectrum made 
from these ^’-factors with a I-A spectral resolution is presented in Fig. 


1. The R { + R Q 1X branch occurs at 3238 A, consisting of about 15 
individual lines. In Fig. 2 we present ^-factors for this branch as a 
function of heliocentric velocity al r - 1 AU. The ^-factor for this 
branch is 6.2 x 10’ 6 sec" 1 for u(r) = 0 km sec - 1 and 7.8 x 10 ft sec -1 
for u(r) = - 28.5 km sec" 1 . The ^-factor of the entire 0-0 band ranges 
from 2.5 x 1 0 _ 5 to 4 . 1 x 10“ 5 sec" 1 for the heliocentric velocity between 
-80 and 80 km sec -1 at a constant temperature of 200 K. The total 
factor also varies between 3.4 x 10" 5 and 4.4 x 1 0 - 5 see -1 for the 
temperature range of 50-400 Kat a heliocentric velocity of 0.0 km sec -1 . 
Our values are significantly lower than the value (2.32 x 10“ * sec -1 ) 
derived by Krishna Swamy and Wallis (1988), who apparently did not 
consider predissociation in the A-X band. 

Krishna Swamy and Wallis (1986) and Wallis and Krishna Swamy 
(1987) tentatively identified the B-X band of SH in the 1670-1920 A 
range in the IUE spectra of certain comets. Krishna Swamy and Wallis 
(1988) derived the ^-factors of the B-X bands to be less than 5 x 10’ 8 
sec -1 , which is so low that the presence of the band in the noisy IUE 
spectra is doubtful. They also did not consider possible predissociation 
in the B-X band. According to the absorption spectrum of the B-X band 
presented in Fig. 1 of Morrow (1966), the absorption linewidths (-0.8 
cm -1 ) are wider than the uncertainty (±0. 1 cm -1 ) in the line positions. 
This suggests that predissociation occurs through these lines, and the 
predissociation will lower the derived ^-factors even more, rendering 
the tentative detection even more doubtful. 

Upper limits of SH in Comet P/Brorsen-Metcalf and lifetimes of 
SH. In order to calculate the upper limits of SH in comets, the lifetimes 
of H 2 S and SH are required. Crovisier et al. (1991) estimated an H 2 S 
lifetime of 4000 sec at 1 AU using an absorption spectrum of Lee et al. 
(1987). 

The dissociation rate of SH through the 0-0 transitions can be easily 
derived from Eq. (1) and is given by 

£> r (0-0) = 'Zx ar B rfPrr . (4) 

/" 

Here we have used the lifetimes given previously to assume that the 
branching ratio for dissociation is nearly unity. The total dissociation 
rate through the 0-0, 1-0, and 2-0 bands is the sum of Z) r (0-0), D r (1-0), 
and D r (2-0). D r (0-0) at 1 AU is 8.20 x 10" 3 sec" 1 for v(r) = -28.5km 
sec -1 and 5.78 x 10“ 3 sec -1 for v(r) = 0 km sec’ 1 . According to 
Franck-Condon (F-C) factor calculations by Nicholls et al. (I960), 
which, although old, is the only available source, the F-C factors of the 
0-0, 1-0, and 2-0 transitions are 0.74, 0.20, and 0.04, respectively. The 
0-0, 1-0, and 2-0 bands occur at approximately 3250, 3070, and 2923 
A, where the solar radiation densities are about 1.2 x 10' 3 °,6.0 x 10’ 21 , 
and 5.2 x I0" 21 erg cm -3 Hz -1 , respectively. Considering the F-C 
factors and the radiation densities, w ; e obtain a total dissociation rate of 
9.50 x 10“ 3 sec -1 at I AU for u(r) - —28.5 km sec’ 1 . The SH lifetime 
is, therefore, about 105 sec at 1 AU for v(r) = -28.5 km sec’ 1 . 

O'Dell etal. ( 1991 ) presented a spectrum of Comet P/Brorsen-Metcalf 
covering the 3000-3600 A range with 1-A spectral resolution. From Fig. 
4 of O'Dell et al ., we obtain an upper limit flux of 5 x 10“ 17 erg cm’ 2 
sec" 1 pixel" 1 for the 3238- A band of SH. Since this flux is an average 
over the entrance slit, 1.28 by 139 arcsec (1.28 arcsec per pixel), the 
total flux, F, within the slit is 5.4 x 10" 15 erg cm" 2 sec’ 1 . At the time 
of the observations, the average heliocentric distance was 1.0 AU, the 
average geocentric distance was 0.645 AU, and the average heliocentric 
velocity was -28.5 km sec -1 . Using the formula 

MSH per slit) < — f~ f (5) 

ghv 

we obtain N (SH per slit) < 1.3 x 10 29 for the slit area of 600 x 32,600 
km. Here, A is the geocentric distance, j? is the ^-factor, and v is the 
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FIG. 1. A model spectrum of the ()-() band of the A- X system with a 1-A slit width for a heliocentric velocity, u(/*) = -28.5 km sec at 
r = I AU. 


frequency. Using the H 2 S lifetime of 4(K)0 sec and the SH lifetime of 105 
sec, and a Haser model (Haser 1957) with u = I km sec _! for SH and 
H : S, we obtain a production rate of SH. Q (SH), which is less than 
6.7 x 10 27 sec 1 . Using the OH (0-0) observations of Comet B-M by 
O'Dell et al., we determined Q (H : 0) to be about 4.0 x 10** sec " 1 using 
a Haser model with v = 1 km sec 1 . We derive that Q (SH)/(9 (H 2 0) is 
about <0.017. 

Crovisier et al. (1991) reported a production rate of H : S, <2(H : S)/ 
£)(H 2 0) = 2 x l()\ for both Comets Austin and Levy (1990c) using 
their microwave observations. Our upper limit for Comet B-M is thus 
consistent with the results of Crovisier et al. for Comets Austin and 
Levy. Detection of SH, which should require an order of magnitude 


greater sensitivity, would be valuable in testing the models used to 
interpret the observations of H : S. 
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